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Abstract 

In this work we propose an approach to deal with radiation field states which incorporates damping effects 
at zero temperature. By using some well known results on dissipation of a cavity field state, obtained 
by standard ab-initio methods, it was possible to infer through a phenomenological way the explicit form 
for the evolution of the state vector for the whole system: the cavity-field plus reservoir. This proposal 
turns out to be of extreme convenience to account for the influence of the reservoir over the cavity field. 
To illustrate the universal applicability of our approach we consider the attenuation effects on cavity-field 
states engineering. The main concern of the present phenomenological approach consists in furnishing a 
straightforward technique to estimate the fidelity resulting from processes in cavity QED phenomena. A 
^ ' proposal to maximize the fidelity of the process is presented. 
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After more than half a century since their proposition by Einstein, Podolsky, and Rosen (EPR) Q, entangled 
states have become the cornerstone of a set of striking proposals in theoretical physics. The advent of Bell's theorem 
Q , which has permitted to test empirically the phenomenon of nonlocality, and more than a decade of experiments 
confirming such an astonishing character of quantum mechanics gave to entanglement the necessary credibility to 
motivate such proposals. Ranging from quantum communication 0] and cryptography to teleportation |^|,^| and 
quantum computation |Q , the nonlocal character of entangled states has let its original realm of quantum metaphysics 
to inaugurate the possibility of a new technology. In fact, inspired by these theoretical implications, the experimental 
setups, formerly designed to investigate fundamental physics such as quantum coherence and nonlocality, have recently 
been enhanced for the realization of teleportation [|| and to demonstrate quantum-logic operations Q ■ 

However, the above-mentioned processes on quantum mechanics come up against a crucial problem intrinsic to 
quantum nature: the decoherence coming from the inevitable coupling between quantum systems and the surrounding 
environment. Such a decoherence process, transforming a pure state in a statistical mixture p, constitutes the 
main difficulty preventing the realization of the above mentioned quantum processes. Since the transformation |?/>) 
— > p becomes faster as the excitation of the initial quantum state increases, the realization of massive computation, 
for example, turns out to be prohibitive and even the realization of a single logic gate operation || leads to a degraded 
output. In this connection, the calculation of the fidelity of a given process becomes a crucial task when proposing 
its experimental realization. However, this task turns out to be as difficult to achieve as the implementation of the 
protocol for accomplishing a given quantum process. As an example, we cite the process of engineering a cavity-field 
state, which requires a series of steps to be pursued ||, since in most of these schemes the cavity- field is built up 
photon by photon. 

To overcome the difficulty in estimating the fidelity in cavity QED processes, in this paper we propose a phe- 
nomenological way to handle with dissipation of a cavity-field state under the influence of a reservoir at absolute zero 
temperature. Our strategy is to incorporate, in an concise algebraic approach, the main results obtained by standard 
ab-initio techniques on dissipation of a cavity field. In short, we include the exponential decay law accounting for the 
damping of the cavity-field excitation explicitly in the evolution of the whole system comprehending the cavity field 
plus reservoir. Once we have calculated the evolved whole state vector, it is immediate to obtain the total density 
matrix, and by tracing out the reservoir variables we finally get the desired reduced density matrix for the cavity 
field, which is the same as the one obtained by standard methods. This approach resembles the Monte Carlo Wave 
Function method pXj ] in the sense that we work directly with the wave function, providing an efficient computational 
tool to account for dissipation in quantum optics. 

It is worth mentioning that our result for the evolved state vector of the radiation field plus reservoir can be applied 
to calculate dissipation effects in any physical process based on cavity QED. Here, as an application, we consider 
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damping effects in engineering a quantum field state and present an alternative method to maximize the fidelity of the 
engineered state. The process of generation of an arbitrary radiation field state, as will be here analyzed, is severely 
damaged due to the reservoir attenuation. 

This paper is organized as follows. In section II we describe the essential features of the model here considered, a 
monochromatic cavity-field trapped in a lossy cavity. In section III we present an algebraic way to introduce damping 
explicitly on the evolution of the whole state vector for the cavity field plus reservoir. In section IV we apply the 
technique developed in section III to analyze the fidelity of an engineered cavity-field state. Section V is devoted to 
illustrate a technique to optimize the fidelity of the engineered field state and finally, in section VI we present our 
conclusions. 



II. MODEL 



We are concerned with the standard model of a cavity field described by creation and annihilation central oscillator 
amplitude operators a', a, coupled to a reservoir consisting of a set of oscillators described by amplitude operators 
b\, bk- The reservoir simulates the cavity damping mechanism: photon losses either at the mirrors or by leaking out 
of the cavity. The cavity field is coupled to the reservoir through the usual rotating wave approximation, so that the 
Hamiltonian for the system is written as 

H = Titoa^a + ^2 Tiu k b\b k + ^ fi (^gkab\ + gla^b^j , (1) 

k k 

where g^ are coupling parameters. The Heisenberg equations for both amplitude operators a and b are given by 

a (t) = -iuia- i^gkbk, (2a) 

k 

bk (t) = -i^kbk - ig* k a. (2b) 

Considering the Wigner-Weisskopf approximation, in which the frequencies of the reservoir oscillators densely cover a 
range in the neighborhood of the field frequency, we obtain the solutions of the equations corresponding to the system 
(^a|,^l]). The solution for the cavity field reads 

a(t)=/*(t)a(0)+53t?fc(f)&k(0). ( 3 ) 

k 

Introducing the damping factor 7 and the typical small frequency shift Alu, the function /i(t) is explicitly given by 

/i(t) = exp j- ^+i(Alu + lu) tj. 

For the purpose at hands it is not necessary to specify the functions $k (t) ■ 

A convenient way to calculate the reduced density matrix for the cavity field consists in using the standard char- 
acteristic function technique ]l2|] . For this intent, we must first calculate the characteristic function \i then the 
Glauber-Sudarshan P-representation, and finally the reduced density operator. The characteristic function can be 
written in terms of the operators a, in the normal order form as 

X (r), rf,t) = Tr {p(t) exp [r,J (0)] exp [-ry*a(0)] } (4a) 
= Tr{p(0)exp [rja^t)] exp [-r)*a{t)\ } . (4b) 

Eq.([4a|) refers to the Schrodinger picture and the Eq. ( fih| ) to the Heisenberg picture. Since the time evolution of 
operators a,a^ are known from Eq.(|J) together with the initial state of the whole system, Eq. (flb|) turns out to be 
more convenient to calculate the characteristic function. Considering that the cavity field is initially in an arbitrary 
superposition state J2n=o^ n wnere \ n ) IS a "-photon Fock state, and the reservoir (assumed at absolute zero 
temperature) is in the vacuum state, the total density matrix at t — is given by 

N 

P (0)= C n C m \n) H <8> |{0}) <{0}| , (5) 

n.m— 
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where we have denoted Y\ k |0fc) = |{0}). 

The Glauber-Sudarshan P-representation of the density matrix can be written as a two-dimensional Fourier trans- 
form of the characteristic function \, 

P(a,a*,t) = J ^x(?7,ry*,i)exp(-77a* +rj*a) , (6) 

while the reduced density operator for the cavity field is given by, in the general form, 

p(t) = I d 2 aP(a, a*,t) \a) (a\ . (7) 



Considering all these definitions it is possible to obtain, after a little algebra, the explicit form for the reduced 
density operator of the cavity field 



l\ (m — j)\ 

n,m=0j=0 1=0 v J ' 



Im.l 

\j + n-l-m)(j-l\. (8) 



(J + n-l-m)\(j-l)\ 



III. PHENOMENOLOGICAL APPROACH 



The aim of this section is to include explicitly the energy dissipation of the cavity field in the time evolution of the 
state vector for the whole system. In short, the exponential decay law accounting for the damping of the cavity-field 
excitation will be explicitly included into the evolution of the state vector for the whole system. Such a state vector 
must lead, naturally, to the result derived from formal techniques for the reduced density matrix Eq.(|^). As usual, 
the reduced density matrix for the cavity field is obtained by tracing over the reservoir variables, i.e., 

ps(t) = Tr R [ps+R(t)} , (9) 

where ps+H.(t) is given, for pure states at time t, by 

PS+R(t) = \*s+R(t))(*s+R(t)\, (10) 

and the subscripts S and R refer to the system and reservoir, respectively. What is thus necessary to be known 
is precisely the evolution of the whole state vector |\I/s_|_ij (t)) . For this intent, we begin by considering first some 
examples which will guide us to the form of the evolution of the generalized n— photons cavity-field state in the 
presence of the reservoir. For the simplest case, a single photon in the cavity, the evolution of the whole state vector 
can be phcnomcnologically inferred as 

|1) \R) — > cxp (-It) |1) \R) + |0) ]T a k b\ \R) , (11) 

k 

where each ctk represents the probability amplitude that a photon have been absorbed by the fc-th mode of the 
reservoir, and 7 is the damping rate for the one-photon state in the cavity. Under the normalization condition on 
state ([ll]), we arrive at the constraint 

^K| 2 = l-exp(- 7 ;). (12) 

fe 

For the next case of a two-photon state, we infer the evolution of the whole system by noting that the interaction 
between the cavity mode and the reservoir is always characterized by one-photon exchange. In this way, the final 
form of the evolution of the whole state vector can be guessed as 

|2) \R) — cxp (-2^) |2) \R) + 

V2exp (-^*)|l>5>* 6 fcl fl > + 

k 

\Q)Y,<* k a k >blbl\R). (13) 

kk' 
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The first term on the r.h.s of Eq.(13) is associated with the probability that the cavity field remains in its original 
two-photon state. It is worth noting the difference between the relaxation rates for one- and two-photon field states: 
the damping rate for a two-photon state is increased by a factor 2 compared to that of the one-photon state. In the 
second term on the r.h.s of Eq.(|l3|), we have included the probability that the one-photon state has originated from 
the decay of the two-photon state and also the probability that the one-photon state remains as such. The last term 
accounts for the probability that all photons have been absorbed by the reservoir. We stress that the imposition of 
the normalization condition on state Eq . (|l3|) gives the same constraint Eq. (|l2]) . 

We note that states Eqs.(|Tl]) and Jl3|), as it should be, recover the same result as given by Eq.(||). To verify this 
assertion, it is necessary to build up the density operator for the system plus reservoir and after that to get rid of the 
reservoir variables. 

Finally it is possible, by induction, to generalize the above results to construct the evolution of an arbitrary n-photon 
state. This state vector for the whole system includes the ingredients discussed before and can be read as 



\n) \R) — > exp (-n~i) \n) \R) 



l!(n-l)! 



exp 



(n- l)M 
V ; 2 



■ exp 



2!(n- 2)! 
+ |0) c kl ...a kn bl...bljR) 



(n - 2) It] \n-2)J2 ^ k2 b\b\ 2 \R) 



(14) 



where each term appearing in Eq.(]l4|) includes the probability that a given state has originated from the initial n- 
photon state and also the probability that it will remain the same. Again, the desired reduced density matrix for the 
cavity field can be obtained from the time-dependent state vector Eq.(|l4|) by using Eqs.(^) and (pX|). 



IV. NOISE EFFECTS ON ENGINEERING QUANTUM-FIELD STATES 

To illustrate the applicability of the method developed in the previous section, we restrict our analysis to the 
damping effects in the engineering process of a cavity- field state. To do that, we choose Vogel et al.'s scheme to 
engineer an arbitrary cavity-field state [Tl| ] , which is based on successive resonant interactions of M two- level atoms 
with an initially empty cavity. The experimental setup for this scheme is depicted in Fig.l. The Rydberg atoms are 
laser excited before entering into the Ramsey zone R, placed in the way of the atoms to the cavity C, preparing each 
atom in a particular superposition of excited |e) and ground \g) states, as is required to properly build up the cavity 
field, photon by photon. Each two-level atom , after interacting with a monochromatic field state in the cavity C, 
leaves its photon in the cavity, being necessarily detected in its ground state by detector chambers D (D e and D g for 
ionizing the states |e) and \g), respectively). 

To exemplify an application of the present phenomenological approach, we consider the first step of this process. 
Initially, the cavity field C is in the vacuum state while the atom is in an arbitrary (unnormalized) superposition |e) 
+ie\ \g), where t\ is a complex parameter controlled by the Ramsey zone R. Once we are considering the reservoir 
at absolute zero temperature, it will always be in the vacuum state \R) = |{0}). 

To simplify the calculations, we assume a stable excited atomic state and just take into account the errors introduced 
by the cavity dissipation mechanism. This is a good approximation, since on average eight out of ten atoms are able 
to travel the distance of the whole setup without decaying. In fact, since a Rydberg-atom excited state has a lifetime 
(l/7 a ) of the order of 10~ 2 s [jl3), the probability of staying in this state is about 0.8 for an experiment duration of 
about 2 x 10~ 3 s. For high-Q superconducting cavities, the lifetime (I/7) is also of the order of 10~ 2 s [ fI3| . Finally, 
we neglect the field dissipation during the time that the atoms interact with the cavity field. 

After the resonant atom-field interaction, when the atom has left the cavity and measured in the ground state, the 
state of the whole system can be written as 

1 

where we have described the damping effects on the cavity field (given by Eq.(|ll])) through the operators 
appearing above, which are explicitly given by 
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P^oghl (16a) 



k 



?p = exp 



H*')^' (16b) 



where t' is the time in which the relaxation process occurs, i.e., the time interval in which the first atom left the cavity 
and is detected. At this point it is useful to define the parameters Cn = cos(gTk\/n + 1) and = sm(gTk\/n + 1), 
where is the interaction time of the fc-th atom with the cavity field and g is the atom-field coupling constant, 
assumed to be the same for all the atoms. The constants ip[ , i = 1,2 in Eqs.( 16a| . 16b ) are the same as those in Ref. 
0, i-e., 

<pP=$\ (17) 
<pM = -e x . (18) 

The next step consists in the analysis of the results after the passage of the second atom through the cavity and the 
subsequent relaxation of the radiation field then constructed. After the second atom (prepared in the superposition 
|e) +ie2 \g}) has left the cavity and again been measured in the ground state, the two-photon cavity field-reservoir 
state becomes 

* (2) )=EQl 2) N|i?>, (19) 

71 = 

where the operators Qn^ turns out to be more complicated then those in Eq.([l5|) due to entanglement between the 
states of each atom and the cavity field-reservoir states. These operators can be written as 



k k\k2 

Q? = exp (-It) + V2T f exp (-If) £ af% (20b) 

k 

Qf =f 2 (2) exp(-2|t), (20c) 

where the operators Tn appearing above are given by 

?o (2) =~e 2 Q ( o\ (21a) 
TW=SPQP-e 2 cWQP, (21b) 
f 2 (2) =s[ 2) Qf\ (21c) 

Here t is the time starting after the second atom left the cavity. For simplicity, we have assumed the time interval for 
the detections of the first and second atoms as equal, and that the second atom enters into the cavity immediately 
after the first atom being detected. Although here we are just concerned with the desired \ipd) = J2n = o<i n \n) 

cavity field state, the method described above can be generalized for generating an arbitrary cavity field state ( 
IV'd) = Sri=o dn \n)) in a straightforward (if rather tedious) manner. Naturally, all the informations we need are 
contained in l^ 2 )) according to Eq. (|l9|). For instance, the reduced density operator can be obtained readily by 
using Eq. (|l^ ) and definitions Eqs. (g) and ((!(]) as 

p s (t) =J\f[a |2) (2| + b |1) (1| + c |0) (0| + (22) 
+ (d|2)(l| + /|2) (0|+ ff |l) (0|+/i.c.)], 

where the coefficients appearing above are given in appendix A, and Af = 1 /(a + b + c) is the normalization constant. 

From this reduced density operator, Eq. (p2|), which includes the inevitable coupling between the system and the 
reservoir, it is possible to calculate, in principle, any physical quantity (observable) of interest. 

Ideal case. When damping effects are neglected, i.e., letting the damping constant 7 = 0, Eq. ( p2| ) recovers the 

ideal case, and ps can be written in terms of a pure state |\& s) (*s|, where y&s) — J2 n =o ^™ \ n ) an( l * ne coefficients 



ip„ are the same as those in Ref. O], as can be checked from Eqs. (16a) to fl21c|) 
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V. OPTIMIZING THE FIDELITY OF AN ENGINEERED CAVITY-FIELD STATE 



At this point this formulation is completely general, and can be applied to generate any arbitrary cavity field state. 
However, we stress that there is a crucial difference between the procedure described here and the one in Ref. fLl}] . In 
fact, in the treatment proposed by Vogel et al., it is necessary to solve a polynomial equations for ei and 62, which arise 
from the consistence of the recurrence equations relating the amplitudes <fk ' On the other hand, in the treatment 

described here, in virtue of the reservoir, the quantities Qn in Eq. ( jl9| ) are now operators, instead of simply being 
parameters. In this connection, there is no immediate way to extract useful information from the recurrence relation 
associated between the operators involved. Thus, a first question arising is, how to turn round this problem? We 
propose the following way to solve this problem: consider a specific desired cavity-field state to be engineered and 
calculate the fidelity of this ideal state (free from the reservoir effects) with respect to the state |^( 2 )) in Eq. jis|). 
By doing this procedure, the only unknown variables are t\ and £2, since all the others can be fixed. The trick here is 
to find numerically the values for eiand £2 which maximize the resulting expression for the fidelity. This process can 
be even more optimized by an adequate choice of the interaction parameters gTk. 

Another important ingredient in this discussion is concerned with the probability for successfully achieve the 
engineering process. In fact, the success of the engineering process depends on detecting all the required atoms in 
their ground states. The expression for the probability to detect each atom in the ground state can be obtained in an 
usual way. The distinguished feature here is given by the presence of the reservoir. To do that, we have to consider 
the internal atomic states as well as the constructed radiation-field states being taken normalized. We thus have 
to normalize the states given in Eqs. (Ill© by following the recipe) <I<( fe )^ — > j\f k |vt/( fc )^ where the normalization 
constants read 



(23a) 




(23b) 



The total probability for successfully engineering the desired state reads V = Ylu—x"Pk, where Vk — 1/ (Afc) 2 ac- 
counts for the probability to detect the fc-th atom in the ground state \g). 

To exemplify the whole procedure, we restrict our analysis to the generation of the following radiation-field truncated 
phase-state, 

|* rf > = ^(|0> + |1> + |2». 

The fidelity of the density matrix p s (t) relative to Y&d) is given by 

T ={^ d \ Ps (t)\^ d ), (24) 

where p s (t) is the normalized density matrix. The resulting expression is a little bit involved and must be numerically 
maximized regarding the variables e\ and 62- To optimize the engineering process (p4|), we consider different interaction 
parameters gr k to each atom, which depend on experimental capabilities. In the present realistic scheme for engineering 
a cavity-field state two important features arise: the fidelity T of the engineered state and the probability V for 
successfully achieve the process. Fig.2(a,b,c) show the histogram of the probability V, the fidelity T , and the rate 
1Z = VJ-, which serves as a cost-benefit estimate to choose the best parameters, in terms of the parameters gT\ and 
gT2- For a definite value of the interaction parameters gTk, it is calculated the correspondent values ei and £2 which 
maximize the fidelity. From this figure we note that a larger probability does not necessarily imply a better fidelity. 
This is an important point since the convenient probability V and the fidelity T to be chosen depend on the engineer 
necessities. In Table I some values of gTk and the solutions for t\ and 62 as well as the correspondent quantities 
V, T , and 1Z are exhibited. Fig.3(a-f) show both the histogram of the elements of the reduced density matrix and 
the associated Wigner functions. Fig.3(a,b) show the ideal case (7 = 0); Fig.3(c,d) show the realistic engineered 
state by using the parameters e± and €2 obtained from Vogel et. al scheme taking into account the reservoir; and 
Fig.3(e,f) show the engineered state obtained by our maximization procedure. As can be viewed in Fig.3(e,f), this 
procedure improve the quantum nature of the engineered state, which is revealed by the negative portions of the 
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Wigner function. To plot the figures 2 and 3 we used the realistic experimental parameters 7 = 10 2 s~ 1 |b|] for the 
decay rate and the estimated time interval t = t' = IO/7 [[TiJ for the first and second atoms, respectively, to reache 
the detection chambers after leaving the cavity. 



In this paper we have proposed an alternative way to treat damping effects at zero temperature in cavity QED 
processes. Our phenomenological approach consists in incorporate the main results obtained by standard ab-initio 
techniques on dissipation of a cavity field directly in the evolution of the state vector of the whole system: cavity field 
plus reservoir. In summary, we have included the exponential decay law accounting for the damping of the cavity-field 
excitation explicitly in the evolution of the whole state vector of the system. 

By considering a standard model of a given mode of the radiation field coupled to a collection of iV-modes rep- 
resenting the reservoir, it was possible to infer the time evolution for the whole state vector in a straightforward 
manner. We stress that this result is rather general, and in principle can be applied to whichever quantum process 
as, for instance, quantum communication, logic operations, teleportation, and cavity-field state engineering, among 
others. 

The phenomenological-operator approach to dissipation in cavity quantum electrodynamics here presented consid- 
erably simplifies the introduction of the inevitable errors due to the environmental degrees of freedom when describing 
processes involving atom- field interactions. The development of this technique became possible due to the previous 
formal work in quantum dissipation, from which we have recovered the main results concerning energy loss of a 
trapped radiation field. So, such convenient approach precludes the necessity of performing the usually extensive ab 
initio calculations as the standard master equation, the characteristic functions, or even the path integral formalism. 
It is worth noting that the technique here developed at absolute zero can be extended for a thermal reservoir. 

As an application of the present technique, we have considered the process of engineering a cavity-field state in a 
lossy cavity. In order to estimate how far the generated state deviates from the idealized, due to attenuation, we have 
analyzed the fidelity of the engineering process. Also, we have proposed an alternative way for engineering a cavity 
field state in the presence of a reservoir, which consists in maximizing the fidelity of the desired (ideal) state relative 
to the (damped) engineered state. A specific example has been given which takes into account realistic values for 
the parameters involved. To conclude, we should stress that damping effects seriously restrict such an engineering 
process. However, these noise effects can be minimized by the scheme proposed in this paper. 
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Appendix A 

In this appendix we show explicitly the coefficients appearing in Eq. (E3), 
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"hi 2 (^f)V^(l- e ^) + 

+ (W)'(l-e-^) (l-e^) + 

-| £l | 2 (5f) 2 (l-e-^) + 
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Figure Caption 

FIG. 1. Sketch of the experimental setup for engineering a cavity-field state. 

FIG. 2(a,b,c). Histogram of the probability P, fidelity T, and rate 1Z, respectively, in the griX gr 2 plane. 

FIG. 3. Elements of the reduced density matrix and the correspondents Wigner distribution functions: (a,b) for the 
ideal state ; (c,d) for the realistic engineered state built with parameters ei and e 2 obtained from Vogel et al. scheme 
including the reservoir; and (e,f) for the optimized engineered state by the maximization procedure. The Wigner 
W (a, a*) distribution functions are plotted as a function of q = Re(a) and p = Im(a). 

Table 

Table 1. Some values of the interaction times gr\ and gr 2l Ramsey zones parameters e\ and e 2 (adjusted in order 
to maximize the fidelity), fidelity T, probability V, and rate 1Z. 

gn gr 2 ei £2 T V Tl 

0.6 3.0 2.7693 -0.1583 0.9253 0.0697 0.0645 

1.3 1.2 -0.8508 + ^0.2874 -0.8838 - iO. 1600 0.8789 0.8831 0.7761 

1.4 2.8 1.2349 + i0.8632 -0.3583 + z0.2427 0.9087 0.3637 0.3305 
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